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Abstract 

We study the quark helicity distributions at large x in perturbative QCD, taking into account 
contributions from the valence Fock states of the nucleon which have nonzero orbital angular mo- 
mentum. These states are necessary to have a nonzero anomalous magnetic moment. We find 
that the quark orbital angular momentum contributes a large logarithm to the negative helicity 
quark distributions in addition to its power behavior, scaling as (1 — x) 5 log 2 (l — x) in the limit of 
x — » 1. Our analysis shows that the ratio of the polarized over unpolarized down quark distribu- 
tions, Ad/d, will still approach 1 in this limit. By comparing with the experimental data, we find 
that this ratio should cross zero at x ~ 0.75. 

PACS numbers: 12.38.Bx, 12.39.St, 13.85.Qk 
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1. Introduction. Power-counting rules for the large-x parton distributions in hadrons 
have been derived more than three decades ago based on perturbative quantum chromody- 
namics (pQCD) combined with a S"-wave quark model of hadrons The basic 

argument is that when the valence quark carries nearly all of the longitudinal momentum of 
the hadron, the relevant QCD configurations in the hadronic wave function become far off- 
shell and can be treated in pQCD. A generic factorization has recently been used to justify 
the power-counting rule by relating the parton distributions at large- x to the quark distri- 



bution amplitudes of hadrons |5l. The power-counting rule has also been generalized to sea 

n q 

marks, gluons, helicity-dependent distributions p, 13], and generalized parton distributions 



The leading pQCD diagrams associated with the leading Fock state of the proton wave 
function predict that the positive helicity (quark spin aligned with the proton spin) quark 
distribution q + (x) scales as (1 — x) 3 , (x = k + /P + is the light-cone momentum fraction of the 
struck quark and is identical to the Bjorken xb in the leading twist approximation), whereas 
the negative helicity (quark spin anti-aligned with the proton spin) quark distribution q~ fx) 
is suppressed by (1 — x) 2 relative to the positive helicity one, scaling as (1 — x) 5 at large x 3]. 
The direct consequence of these power laws for the quark distributions is that the ratio of 
polarized quark distribution Aq(x) = q + (x) —q~(x) over the unpolarized quark distribution 
q(x) = q + (x) +q~(x) approaches 1 in the limit x — > 1; i.e., at large x, q + dominates over q~. 
When this prediction is compared to the experimental data 0, [lo|, H, 12]. it is interesting 
to observe that, for the up quark the ratio increases with x, and seems to approach 1 at 
large x. However, the ratio for the down quark is still far below 1, and remains negative for 
a wide range of x < 0.6 [9jj. This discrepancy has stimulated much theoretical interest. 

In this paper we will reexamine the large- x quark helicity distributions in the perturbative 
QCD framework {jj 4]]. We work in light-cone gauge with A + _ = 0, where there is no ghost 



contributions, and orbital angular momentum is physical [13J]. We will take into account 
the contributions from not only the leading light-cone Fock state expansion of the nucleon 
wave function with zero quark orbital angular momentum (L z = 0), but also the valence 
Fock states with nonzero quark orbital angular momentum (L z 7^ 0). These contributions are 
naturally required to obtain a nonzero anomalous magnetic moment for nucleons [14j and are 
also present in the wave function solutions in the AdS/CFT correspondence approach 15]. 



We find that for the negative quark helicity distribution q , there exist large logarithmic 
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enhancements from the \L Z \ = 1 Fock state component of the proton. With this large 
logarithmic modification, we can explain the discrepancy between the power-counting rule 
and experimental data. 

2. Analysis of the large- x behavior of the quark helicity distributions. In the 
x — > 1 regime where the struck quark has nearly all of the light-cone momentum of its parent 
hadron, the relevant QCD dynamics becomes far-off the mass shell: the Feynman virtuality 
of the struck quark becomes highly space-like: k 2 F —m 2 ~ _ ^L±^j ; where k± is the transverse 
momentum of the struck quark and A4 is the invariant mass of the spectator system. Thus 
we can use perturbative QCD to analyze the large-x behavior of the parton distributions 
since the internal propagators in the relevant Feynman diagrams scale as 1/(1 — x). This 
behavior leads to the power-counting rules. In fact, more partons in the hadron's wave 
function mean more propagators in the scattering amplitudes and more suppression for the 
contribution to the parton distributions. Thus the parton distributions at large-x depend on 
the number of spectator partons in the Fock state wave function of the hadron. For example, 
the valence quark distributions of nucleon will be dominated by the three-quark Fock states 
of the nucleon wave function. The three-quark Fock state expansion of the nucleon wave 
function consists of zero orbital angular momentum component (L z = 0) and nonzero orbital 
angular momentum component (L z 0) [16]. In the following discussion, we will consider 
the contributions from both Fock state components. 

In the original power-counting analysis of the quark helicity distributions [6], only the 
contributions from the leading Fock state with L z = have been taken into account. In 
Fig. l(a,b) we show the typical diagrams which contribute to the positive (a) and negative 
(b) quark helicity distributions at large x. In terms of the leading order quark distribution 
amplitude which corresponds to the L z = Fock state expansion of the proton wave function 
{4], these quark helicity distributions can be written as, 

s^Ui = J [dy.iWMyiWiy^ fa, 16 (i - *)) , (i) 

where the integration measure [dz/i] is defined as [dyi\ = dyidy 2 dy 3 5(l — 2/1 — 2/2 — 2/3), and 
the yi are the light-cone momentum fractions of the proton carried by the quarks in the 
light-front wave functions, i.e., Pi = yiP and p\ = y[P in Fig. 1. Here, <3? and $' represent 
the quark distribution amplitudes of the proton at the left and right sides of the cut line, 
respectively. H represents the hard part of the amplitude which depends on y t and y' { , and 
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FIG. 1: Examples of Feynman diagrams which contribute to the quark distributions at large 
x: (a) for q + with power contribution of (1 — x) 3 ; (b) for q~ with (1 — x) 5 ; (c) for q~ with 
(1 — x) 5 log 2 (l — x). The wave functions at the left and right sides of the cut line in (a) and 
(b) represent the leading Fock state expansion with zero quark orbital angular momentum, whereas 
those for (c) represent the valence Fock state with one-unit of quark orbital angular momentum. 

(1 — x) as well. 

The large-x behavior of the hard factor 7i can be evaluated from the partonic scattering 
amplitudes as shown by the Feynman diagrams in Fig. 1. The wave functions corresponding 
to Figs. l(a,b) have zero quark orbital angular momentum, and thus the total quark spin 
will be equal to the proton spin. If the struck quark spin is the same as the proton spin (for 
the positive helicity quark distribution q + ), the two spectator quarks will be in the spin-0 
configuration, whereas for the negative helicity quark distribution q~ with the quark spin 
opposite to the proton spin, the two spectator quarks will be in the spin-1 configuration. It 
has been known for a long time that the hard partonic part TC with a spin-1 configuration for 
;he two spectators will be suppressed by (1 — x) 2 relative to that with spin-0 configuration 
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171 ]. That is also the reason why the negative helicity quark distribution is suppressed 
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by (1 — x) 2 relative to the positive helicity quark distribution from this contribution. 

The far-off-shell propagators in the partonic Feynman diagrams are the dominant contri- 
butions to the power-counting of (1 — x) at large x. One of the gluon propagators in Fig. 1 
behaves as 

1 1 1 l-x 



(P3 - k 2 ) 2 2p 3 ■ k 2 (k 2 ± ) y 3 
at large x. In the above expression, we have omitted all higher order terms suppressed by 
(1 — x). Here, (k±) ~ (kf^) ~ (fcfj_), represents a typical transverse momentum scale for the 
spectator system. Each propagator will provide a suppression factor of — Counting the 
hard propagators in Fig. 1(a), we find the total suppression factor from the hard propagators 
is 

(1 ^ )8 (3) 



2/22/3(1 - 2/2)2/ 2 2/ 3 ( 1 - V2) 
We notice that the above expression does not introduce additional dependence on (1 — x) 

upon integration over y^ and y^, assuming that the leading twist distribution amplitudes <E> 

and $' are proportional to 2/1J/2J/3 and 2/12/22/3 Eh respectively. Combining these results with 

the power behavior for the other parts of the partonic scattering amplitudes and the phase 

space integral, we find the positive helicity quark distribution o + scales as (1 — x) 3 , whereas 

the negative helicity quark distribution q~ scales as (1 — x) 5 

In the above analysis, we only considered the contributions from the leading Fock state 
of the proton with zero quark orbital angular momentum. In general, the contributions 
from the higher Fock states and the valence Fock states with nonzero quark orbital angular 
momentum will introduce additional suppression in (1-x) j^Q]. However, the nonzero-quark- 
orbital-angular-momentum Fock state can provide large logarithmic enhancement to the 
helicity flip amplitudes. For example, it was found that the nonzero quark orbital angular 
momentum contributes a large logarithmic term to the nucleon's helicity-flip Pauli form 
'actor F 2 (Q 2 ), which leads to the scaling behavior F 2 (Q 2 ) ~ log 2 (Q 2 /A 2 )/Q 6 at Q 2 — > 00 
is[ [lol . This is consistent with recent experimental data from JLab [2(3]. In the following, 
we will study the nonzero quark orbital angular momentum contribution to the q~ quark 
distribution which is also associated with the helicity-flip amplitude. The corresponding 
contributions to the positive quark helicity distribution are always power suppressed Q]. 

In Fig. 1(c), we show an example of a contribution from the L z = 1 Fock state of proton. 
Because the quark orbital angular momentum contributes one unit of the proton spin, we 
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can have difference between the total quark spin and the proton spin. If the two spectator 
quarks are in the spin-0 configuration, this will enhance the power-counting in the hard 
factor 7i. On the other hand, in order to get a nonzero contribution, we have to perform the 
intrinsic transverse momentum expansion for the hard partonic scattering amplitudes {19I ]. 



which will introduce an additional suppression factor in (1 — x) For example, one of the 
contributions from the diagram shown in Fig. 1(c) to the negative helicity quark distribution 
will be proportional to 

q -(x) oc [ (K+^D(pr-^)^ (3 H^,^)^ (3 Hi/:,^)^(^,m;i/:,^) , (4) 



where -ip^ is a light-front wave function amplitude for the L z — 1 Fock state of the proton 
. The intrinsic transverse momentum expansion is essential to get the nonzero contribu- 
tions. Otherwise, the integral over the transverse momenta pi± and p' i± will vanish because 
of the explicit factors p\ + ip\ and p'f + ipf in the above equation. One intrinsic transverse 
momentum expansion comes from the propagator we mentioned above, 

1 _ 1 

(P3 - k 2 ) 2 (V3P -k 2 + p 3 ±) 2 

(3(1 -x) ( P(l-x) 



where (3 is defined as k^ /(l — x)P + , and we have kept the linear dependence on p 3 ± in the 
above expansion. Only this linear term will contribute when integrating over pi±: J k 2 ± ■ 
Ps±(Pi+ip\)ip^ cx (kp +ikn)yaQ>4.(y-\ , y 2 , 2/3), where $ 4 is one of the twist-4 quark distribution 
amplitudes of the proton 19|, [2l|] . From the above expansion, we find that this term will 
introduce additional factor of (1 — x)/y^ in the hard factor. Similarly, because of the p'± —ip'f 
factor in Eq. (4), we have to do the expansion in intrinsic transverse momentum associated 
with the wave function at the right side of the cut line, and again the expansion of the gluon 
propagator with momentum of p' 3 — k 2 will introduce another suppression factor of (l — x)/y' 3 
in the hard factor. Thus the total suppression factor from the above two expansions will 
be (1 — x) 2 /y^y 3 , which gives the same power counting contribution to q~ as that from the 
leading Fock state with L z = in the above. 

We thus find the contributions from L z — 1 Fock state of the proton do not change the 
power counting for the q~ quark distribution at large x. However, the additional factor 
1 / 2/32/3 from the intrinsic transverse momentum expansions will lead to a large logarithm 
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when integrating over yi and y\. This is because, combining the above two factors with all 
other factors from the propagators shown in Eq. ([3]), the total dependence on yi and y[ for 
the hard factor will be 

~ y 2 yi(l - y 2 )y> 2 y>i(l - y> 2 ) ' (6) 
where we have y\ and y' 3 2 in the denominator. On the other hand, we expect the twist- 
4 quark distribution amplitude to have the following behavior at the end point region: 



2/3^4(2/1, 2/2, 2/3) oc yiy 2 yz and ^$4(2/1,2/2,2/3) oc y'lVWz UM- Thus we wil1 have logarithmic 
divergences for the integrations over y 3 and y' 3 , for which we can regularize in terms of 
log(l — x) as indicated in the above propagator expansion. This will lead to a double 
logarithmic contribution log 2 (l — x) in addition to the power term (1 — x) 5 to the q~ quark 
distribution at large x. 

In summary, for the negative quark helicity distribution q~ at large x, the leading Fock 
state with zero quark orbital angular momentum L z = contributes to a power term (1 — 
x) 5 , whereas the valence Fock state with \L Z \ = 1 contributes to a double logarithmically 
enhanced term (1 — x) 5 log 2 (l — x). So, in the limit x — > 1, the q~ distribution will be 
dominated by the contributions from L z — 1 Fock state of the proton, scaling as (1 — 
x) 5 log 2 (l — x). In the intermediate x range, the sub-leading terms can also be important. 



For example in Ref. [6|, the quark helicity distributions were parameterized by the leading 
and sub-leading power terms and fit to the experimental data. This was later updated 
to account for the latest data in Ref. 2^]. Thus, as a first step towards a comprehensive 
phenomenology, we follow the parameterizations for q + and q~ in Ref. [6] by adding the 
newly discovered double logarithms enhanced contributions, 
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} , (7) 

where the additional two parameters C' u and C' d come from the logarithmic modifications to 
the q~ quark distribution at large x, and all other parameters refer to 6|. In the following, we 
will fit to the current experimental data at large x region with the above parameterizations 
for the valence up and down quarks. 
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FIG. 2: Comparison of the quark helicity distributions Eq. ^ with experimental data (left panel), 
and future projections from J Lab (right panel) as functions of x for up (the upper curves) and down 
( the lower curves ) quarks. The circles are for HERMES data lli l. the triangles up for SLA C 11 a] . 



the triangles down for J Lab Hall-A data f^j, the filled squares for CLAS 



id/ , and open squares 
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for 12 GeV upgrade projection for CLAS. The dashed curves are the predictions from ]22J , and the 
solid ones are our fit results (only the large-x (> 0.3) experimental data were used in the fit). The 



symbols in the right panel show combined projections from all three JLab experiments 



3. Phenomenological applications. In order to demonstrate the importance of the 
new scaling behavior for the negative helicity distributions for the valence up and down 
quarks, we analyze the latest experimental data from SLAC, HERMES and Jefferson Lab, 



including Hall A and Hall B data 
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121 ] . We will keep the original fit values for other 



parameters [22| except the two new parameters: C' u and C' d . We only use the experimental 
data in the large-x region, i.e., x > 0.3, where the sea contribution is not significant. We 
perform our fit at a fixed Q 2 = 4 GeV 2 , and all the experimental data are evolved to this 
scale by using the GRSV parameterization [23 j for the polarized and unpolarized quark 
distributions. The evolution introduces some theoretical uncertainties. 
From our fit, we find the following values for C' u and C' d , 



C' u = 0.493 ± 0.249, C' d = 1.592 ± 0.378 , 



(8) 



which are comparable in size to C u = 2.143 ± 0.137 and C d = 1.689 ± 0.227 in Ref. [22]. 
The minimum of the functional x 2 is achieved at x 2 = 11.4 and x 2 /DOF = 11.4/10 = 1.14. 
We further notice that the additional two terms in Eq. ([7]) do not change significantly the 
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sum rules for the up and down quarks, such as the Bjorken and momentum sum rule, which 
are essential for constraining the parameters in Refs. 0, [^j. For example, they contribute 
~ 4% to the momentum sum rule coming from the quarks. 

In the left panel of Fig. [21 we show the above fit, where we plot the ratio of the polarized 
quark distribution Aq over the unpolarized quark distributions q as functions of x for both up 
and down quarks, compared with the experimental data. From these comparisons, we find 
that the ratio for the up quark Au/u can still be described by the parameterization based on 
the original power counting rule for u + and u~ . This can also be seen from the small value 
of C' u in our fit Eq. (JHJ), with big error bar though. However, for the down quark we have 
to take into account a large contribution from the newly discovered term for the negative 
helicity distribution d~; the difference between our result and the original power-counting- 
rule inspired parameterization 22] becomes significant at x > 0.5. The analysis of the 



anomalous magnetic moment and generalized parton distributions of nucleons also indicates 
significant contributions from the orbital angular momenta of up and down quarks 25]. 
This is in qualitative agreement with our fitting results, taking into account the large error 
bar for C' u . A precision determination of these contributions shall be obtained by further 
development for a consistent set of parameters for Eq. (jTJ) from next-to-leading-order QCD 
analysis of both polarized and unpolarized data over the full range in x |22j ]. 

Another important prediction of our fit is that the ratio of Ad/d will approach 1 at 
extremely large x, and it will cross zero at x ~ 0.75. It will be interesting to check this pre- 
diction in future experiments, such as the 12 GeV upgrade of Jefferson Lab. For comparison, 
in the right panel of Fig. 2, we show the experimental projections for these measurements 
from the 12 GeV upgrade of JLab experiments 24|, together with our predictions and the 
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results from the previous power-counting-rule parameterizations 
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